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ABSTRACT 

A macronova (or kilonova) was observed as an infrared excess several days after short gamma-ray 
burst, GRB 130603B. Although the r-process radioactivity is widely discussed as an energy source, 
it requires huge mass of ejecta from a neutron star (NS) binary merger. We propose a new model 
that the X-ray excess gives rise to the simultaneously observed infrared excess via thermal re-emission 
and explore what constraints this would place on the mass and velocity of the ejecta. This X-ray- 
powered model explains both the X-ray and infrared excesses with a single energy source by the central 
engine like a black hole, and allows for broader parameter region, in particular smaller ejecta mass 
~ 10” 3 — 10 -2 Mq with iron mixed as suggested by general relativistic simulations for typical NS-NS 
mergers, than the previous models. We also discuss the other macronova candidates in GRB 060614 
and GRB 080503, and implications for the search of electromagnetic counterparts to gravitational 
waves. 

Subject headings: - 


1. INTRODUCTION 


What is the energy source of a macronovg@ (or kilo- 
nova)? Macronovae are considered as the emission from 
the ejecta with mass ~ 10 -4 — 10 _1 Mq and the velocity 
~ 0.1 — 0.3c, whic h are accompanied by the mergers of a 
NS binary FI (e.g., Li fe Paczvilski lll998HKulkarni 1120051: 


Hotokezaka et al. 120 1 34 IKvutoku. Ioka fc Shibata I 

2013: iKvutoku et al. 112015 ). A NS binary merger is one 
of the most promising sources for the direct detection 
of the gravitational wave (GW). The identification 
of the electromagnetic counterpart to a GW source 
(|Metzeer fc Bergerll2012t iRosswoglI2015D would signif¬ 
icantly reduce the huge localizatio n error o f th e GW 
detectors (~ 10 — 100 deg 2 ; e.g., iBerrv et, al. 1 120151: 
Essick et al. Il201~5h such as Advanced LIGO (Aasi et akl 
20151) , A dvanced VIRGO (lAcernese et al. I I2015D arid 
KAGRA (lAso et al. 1 120131) . The strategy of follow-up 
observations of GW sources should be improved by 
clarifying the main energy source of macronovae, which 
determines the behavior of the light curve. 

Nuclear heating due to the decay of the r-process 
elements has been widely di scussed as a heating 
of macron ovae (e.g.. iLi fc Paczvriski I 119981: 
1201a IKasen. 


source 


Metzger et al. 


Bad nell fc Barnes 
20131 ITanaka fc Hotokezaka I 120131: iRosswog et al 
2014 : ITanaka et al. I 12014 iWanaio et, al. I 12014 
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iLippuner fc Roberts I 12015 ). Since merger ejecta 
are neutron-rich, r -process elements could be synthe¬ 
sized, and the NS binary merge rs co uld be t h e or igin 
of r-process elements (e.g., lLattimer fe Schramm I 
H97l . although iron-rich ejecta are also possible in 

some cases via_the_ sho_c_k and/o r neutrino heating 

(e.g., iMetzger. Piro fc Quataertl 200! );. The r- 


■rocess elements also affec t the o pacity o f t he ejecta 

Kasen, Badnell & Barnes 20131: iThnaka fc Hotokezaka I 

I2013D . making the long timescale and low temperature 
of macronovae. 

The activities of the central engine are 
also proposed as a lternative energy sou rces of 
t he macr on ovae (fKisaka. Ioka fe Takarni I 120151 
iKisaka. Ioka fe Nakamura I I2015D . After a NS binary 
merger, either a BH or a NS is formed as a remnant. The 


ativistic jet (e.g., Kisaka & Ioka 2015), disk wind (e.g., 

Proga & Zhang 

2006^ Kiuchi. Kvutoku & Shibata 

2012; Kiuchi et al. 

[2014L 2015) and NS wind (e.g., 

Dai et al. 20061 

Metzger. Ouataert & Thomoson 

2008fc Rowlinson et al 

20131 Gomnertz et al. 20131 

Gomnertz, O’Brien & Wvnn 2014 Lii et al. 20151). 


which may be observed as prompt, extended and 
platea u emissions in short gamm a-ray bursts (GRBs) 
(e.g., iNakar 1 120071 : iBerger N2014D . These outflows can 
heat the ejecta, which may shine as the macronovae. 
Similar energy inj ection is considered in th e core-collapse 
supernovae (e.g., IKasen fc Bildsten~ll2010l) . The energy 
injection by a highly magnetized NS (magnetar) could 
also produce a br ighter transient than the observed 
macr o novae Ce.g.. lY u Zhan g fc Gao 1 120131: iF an et, al. I 
20131: IMetzger fc Piro I 12014 iWang. Dai fc Yu I 120151 : 


any radio signatures of magnetars (e.g 

.. Nakar & Piran 

2011 

Piran, Nakar & Rosswog 

201! 

Metzger & Bower 

2014 

Takarni. Kvutoku & Ioka 

20141). 


Recently, a macronova candidate is detected as 
the i nfrared excess in GRB 130603B dTanvir et al. I 
120131: IBerger. Fong fc Chornock I I2013D . From obser¬ 
vations, the luminosity, timescale and temperature 
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are ~ 10 41 erg s _1 , ~ 7 days and < 4 x 10 3 K, re¬ 
spectively. In the r-process model, relatively large 
ejecta mass for the merger of a NS-NS binary, 
> 0.03M e , is requi red to reproduce t he observed 
infrared excess (e.g.. IBerger. Fong fc Chor nock I [2 013t 
Piran, Korobkin & Rosswogj 120141 : [Grossman et al. ] 


mW- This is important for the equation of state 


of high-density matter and the progenitor of short 
GRBs (BH-NS merger or NS-NS merger) (e.g., 
iHotokezaka et al.l 2013c ). Too large ejecta mass 
in the r-process m odel may imply the e ngine- 
p ower ed model JKisaka. I oka fc Takami 1 120151 : 

IKisaka. Ioka Sz Nakamura 1 l2015i b~ Still the previ¬ 
ous models using the activities of the extended 
and plateau emissions also re qui re t he l arge ejecta 
mass > 0.02M^ fe isaka, oka fo Takami I 120151 : 

IKisaka. Ioka fc Nakamura I 12015 1. In addition, the 
previous models also require ejecta with high opacity 
~ 10 cm 2 g _1 , which is comparable to the r-process 
case, to reproduce the observed long timescale and low 
temperature. 

GRB 130603B also shows the long-lasting mysteri¬ 
ous X-ray component, which is significantly in excess 
of the ext rapolated power l aws based on the optical 
afterglow (jFong et al. I I2014D . The detected luminos¬ 
ity is ~ 10 42 erg s _1 at the observed time ~ 6 x 10 5 
s, which is longer than the timescale of the plateau 
emissi on dRowlinson et al. 1120131: IGompertz et al. [120131 : 
iGompertz. O’Brien fe Wvnn ll20l4lLu et al. 112015 1. and 


could originate from the central engine such as the 
accretion d isk emission due to the fallback mute- 
rial (e.g.. iRosswogl 120071 : iRossi fc Begelmanl 120091 : 
IKvutoku et ah [Ml). Since the central engine is sur¬ 
rounded by the ejecta, we expect the interaction between 
the emitted X-ray and the surrounding ejecta. Then the 
ejecta heated by the irradiation of X-ray may emit in¬ 
frared photons and reproduce the observed infrared ex¬ 
cess according to the ejecta properties. This is without 
need for any additional energy source such as the radioac¬ 
tive decay of r-process elements. Such an X-ray-powered 
model has a significant advantage that the model uses 
a single energy source to explain two mysterious signals 
- X-ray and infrared - which are observed at the same 
time and with similar luminosities. This is in contrast to 
the r-process and previous engine-powered models which 
require different unrelated sources for the infrared and 
X-ray excesses. 

In this paper, we consider the X-ray-powered macrono¬ 
vae to explain the observed macronova candidate, GRB 
130603B. In Section [2] we describe our model and con¬ 
straints on the ejecta properties. In Section^ we present 
the results for the macronova in GRB 130603B. The 
other macronova candidates in GRB 060614 and GRB 
080503 are discussed in Section [4] We present the dis¬ 
cussion and summary in Section [5] 


2. X-RAY-POWERED MODEL 

Figure CD shows the observed X-ray (1 keV; black), 
optical (r-band; blue) and infra red (H-ban d; red) 
light curves of GRB 130603B dTanvir et al. I 2013 : 
Berger. Fong fe Chornock 1120131: iCucchiara et al. I 2013 : 
de IJgarte Postigo et, al. 112014 iFong et, al. 1 120141 1 . The 
o ptical light curve shows the jet break at ~ 4 x 10 4 s 
(IFong et al. Il2014l j . which is supposed to be achromatic. 



Fig. 1.— Observed light curves of GRB 130603B. Black, blue 
and red denote the X-ray (1 keV), optical (r-band) and in- 


model are taken fromIFon 

2 ; et al. (I2014J) (some of them are orig- 

inated from Tanvir et al. 

20131: Berger. Fong & Chornock 20131: 


the optical and i nfrare d flu xes c orrected by the extinction of the 
host galaxy in lFong et al. 1 620141. The green dot-dashed line shows 
L oc t ~ 2 , on which the X-rays give the same contribution to the 
macronova emission. After the jet break time t > 4 X 10 4 s, the 
X-ray and infrared fluxes are in comparable excesses of the extrap¬ 
olated power laws based on the optical emission. 


However, X-ray and infrared fluxes are significantly in 
excess of the extrapolated power laws after ~ 2 x 10 5 s. 
Remarkably, the X-ray luminosity is similar or slightly 
larger than the infrared excess known as the macronova 
emission at the time ~ 10 6 s. 

The similar luminosities suggest the same origin for 
the X-ray and infrared excesses: the infrared nracronova 
could be reproduced if a part of the X-ray excess is 
converted to the infrared emission, whatever the X-ray 
excess is. This X-ray-powered model is naturally ex¬ 
pected because the NS binary mergers would eject mat¬ 
ter of mass M c j 10~ 3 — 10 as shown by general 
relativistic simulations (e.g.. IHotokezaka et, al. I 2013b; 


Kvutoku. Ioka fc Shibata I 20131 IKvutoku et al. 1201*5 . 

Kawaguchi et al. I I2015D . and provide a screen to ab¬ 


sorb the X-ray and re-emit the infrared emission. This 
model explains the observed infrared macronova with¬ 
out the nuclear heating d ue to the /3-decay and fission of 
r-process elements (e.g., iLi fc Paczvnski 1119981) . 

Note that the activities of the extended emis¬ 
sion (timescale ~ 10 2 s) and plateau emission 

(timescale ~ 10 4 s) may also contr ibute en¬ 

ergy to the macronova ([ Cisaka. Ioka fc Takami 1 120151 : 
IKisaka. Ioka fc Nakamura 120151 ). In these cases, the en¬ 
ergy is injected with shorter timescale than that of the 
macronovae (~ 10 6 s) and remains trapped with adia¬ 
batic cooling (l?i n t oc F -1 ) before being released at ~ 10 6 
s. As a result, the required luminosity for explaining the 
macronova scales as L ~ Flint/l oc t~ 2 . In Figure [Q we 
plot L oc t~ 2 by the green dot-dashed line. Because the 
X-rav excess component follows L oc f -188 (IFong et al. I 
120141) . the later energy injection is more effective for the 
macronova emission. 
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For the X-ray-powered model to produce the X-ray and 
infrared excesses detected by observers simultaneously, 
we consider the following physical setups in Figure [2] 
First, X-rays are generated near the central source and 
are emitted in nearly isotropic direction. Second, the 
ejecta lie at radius larger than the X-ray source and cover 
a fraction of solid angle. Third, the line of sight to the 
X-ray source is clean for the observers who detect the 
GRB emission. 

Although our model does not depend on the specific 
mechanism of X-ray emission, we assume that the X-ray 
excess at ~ 1 — 10 days in GRB 130603B originates from 
the activity near the central engine, such as the accretion 
disk with sup er-Edding ton accre tion rate (e.g.. iRosswog 
2007 : iRossi fc Begelmanl 120091: Fernandez fe Metzger 
20131; iFernandez et al. I l2015aH bl: IKvutoku et al. I 120151) 
like ultra-luminous X-ray sources. At the early stage of 
the merger, a relativistic jet is launc hed from the central 
engine to penetrate the eject a (Nagak ura et al. 1 12014 
Murguia-Bcrt hler et al. If20l4) . Because the jet makes a 
hole in the ejecta, the observers toward the jet axis can 
directly see inside the ejecta (see Figure [2]). 

Let us first consider nearly isotropic ejecta with a con¬ 
stant velocity v e \ base d on the results of the n umeri- 
cal simulations (e.g., Hotokcz aka et alH 1201 3b). The 
anisotropy of the ejecta will be discussed in Section 12.51 
The radius of the ejecta is 


i?n 




(1) 


where t is the time after the merger. Since the velocity 
structure becomes homologous and the ejecta spreads to 
the radial direction, the typical mass density of the ejecta 
is given by 


3 M ej 

fc '~4 rfV 


( 2 ) 


For the composition of the ejecta, we consider both 
the iron-rich ejecta and the heavy r-process ejecta. 
Since the tidally ejected matter is neutron rich, the 
heavy r-process elements may be synthesized (e.g., 
ILattimer fc Schramm 111974 ). On the other hand, the 
iron-rich ejecta may be dominant if the shocks and/or 
neutrino irradiation make the electron fraction high (e g. , 
Metzger. Piro fe Quataertl l2009t iSekiguchi et al. I [20151 : 
Martin et al. ll2015HRichers et al. 1120151) . 

We consider four conditions to reproduce the infrared 
excess in our X-ray-powered model. First, the X-ray pho¬ 
tons should be absorbed by the ejecta, i.e., the optical 
depth for the X-ray absorption rx, a bs should be larger 
than unity, 


7X,abs ^ 1- (3) 

Second, the absorbed X-ray photons should be thermal- 
ized in the ejecta to produce infrared photons. Since the 
opacities for the iron and r -process eleme nts are decreas¬ 
ing f unction of wave length (jKasen. Badnell fc Barnesl 
2013; iTanaka fc Hotokezaka 1 1201311 . we require the op¬ 
tical depth for infrared photons tir to be larger than 
unity, 

tir >1. (4) 





Infrared 

emission 


AAA, 

■ tx. abs ^ ^ 

*IR> 1 
Lliff < t 

'IT < T max 

n. 


Fig. 2.— Schematic picture of the X-ray-powered model. 


tir satisfies condition ([4j) , we consider random walk for 
the propagation of infrared photons in the ejecta. The in¬ 
frared photons can escape from the ejecta if the diffusion 
timescale tdiff is smaller than the dynamical timescale t , 

tdiff < t. (5) 

In the case that the central engine is active longer than 
the photon diffusion timescale, the light curve of the ther¬ 
mal emission follows the evolution of the central engine, 
i.e., the peak time is not necessarily when t ~ tdiff in 
the X-ray-powered model. This is the reason why the 
X-ray-powered model requires smaller amount of ejecta 
than the previous models (see section 3). Finally, the 
temperature of the ejecta T should be lo wer than t he ob¬ 
served upper limit T maY ~ 4000 K (|Tanvir et al. I [20131 : 
iBerger. Fong fc Ghornock~ll2013l) . 

T < T max ~ 4000 K, (6) 

which is obtained from the detected infrared flux and the 
optical upper limit (blue arrow in Figure [TJ). In following 
Sections 12.11 - 12.41 we consider these four conditions to 
give constraints on the ejecta properties. We also con¬ 
sider the luminosity ratio of X-ray and infrared excesses 
in Section 1251 


2.1. X-ray Absorption 

In this subsection, we consider condition Q for the X- 
ray photons to be absorbed in the ejecta. The absorption 
optical depth is given by 

M ej 

T X,abs ~ K bf > ( ‘ ) 

where Kbf is the bound-free opacity. The opacity 
Kbf depends on the ionization state of the ejecta 
(jMetzger fc Piro 112014 . 


«bf 


fn 

Am p 


^"bf 5 


( 8 ) 


Third, the infrared photons should escape from the ejecta 
to be detected by the observers. Since the optical depth 


where A is the average mass number of the composed 
elements of the ejecta, m p is the proton mass, Obf is the 
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photoionization cross section and /„ is the neutral frac¬ 
tion of the ejecta. The ionization state is determined 
by the balance of the photoionization due to the X-ray 
emission from the central engine and the recombination. 
The neutral fraction f n is described by the ratio of the 
absorption rate of the ionizing photons lZ lon and the re¬ 
combination rate lZ re c as 

<•> 

The recombination rate is approximately described by 

R-rec n e a lec , (10) 

where n e is the number density of the electrons and a rec 
is the recombination coefficient. The ionization rate is 


R-ion ^ ph C"bf C- 


( 11 ) 


where the number density of photons whose energy is 
larger than the ionization threshold energy is, 


47r/ir'ii?p i c’ 


( 12 ) 


where hv ; is the ionization threshold energy and e is the 
fraction of the luminosity above the ionization energy. 

In the case of iron-rich ejecta, the number den¬ 
sity of the electrons in the fully-ionized state is 
n e ~ 26p e j/56m p , the recombination coefficient for the 
hydrogen-like iron is a rec ~ 6.3 x 10 -lo (T/10 4 K) -0 ' 8 cm 3 
s _1 (iOsterbrock fc Ferland 1120061) . the ionization thresh¬ 
old energy of the innermost electronic state is hv\ ~ 9.277 
keV, and the photo-ioniz ation cross sectio n at v = v\ i s 
(ibf ~ 1.2 x 10~ 2 ° cm 2 (IOsterbrock fc Ferland 112006 1. 
Using these values and the irradiated X-ray luminos¬ 
ity eLx ~ 10 41 erg s _1 0, the ratio 7^i on /7^rec is much 
smaller than unity for the reasonable parameter ranges 
of the merger ejecta, v e -Jc > 10” 2 , M e j > 10 _3 M e and 
T < 4000 K, at the time t ~ 7 days. Then, we use / n = 1. 
Substituting Equations ©, m, m and (jHJ) with / n = 1 
into condition (131) , we obtain the constraint on the ejecta 
velocity 


— <0.79 

c 




Mej 

10^ 2 Mq 



As long as the neutral fraction f n = 1 is a good approxi¬ 
mation, inequality (fT31) does not depend on the temper¬ 
ature. Note that since the photo-ionization cross section 
for electrons bounded at o uter shell _is la rger t ha n th at 
for the hydrogen-like iron (IOsterbrock fc Ferland 1 12006 ). 
lower energy ionizing photons (< 10 keV) are also ab¬ 
sorbed by the ejecta as long as condition m is satisfied. 

For the r-process ejecta, we assume that most X-ray 
photons are absorbed by the ejecta within ~ 10 days. 


6 We use the X-ray luminosity at 0.3 — 10 keV and re¬ 
gard the excess component as non-thermal emission because the 
photon index is ~ —2 from observational data at ~ 10 5 s 
(http://www.swift.ac.uk/index.php) when the excess component 
dominates the afterglow emission. Since the number of photons 
is smaller in higher energy in the non-thermal emission, we con¬ 
sider that the the photons with higher energy than 9.2 keV do not 
significantly change the ionization fraction. 


If the X-ray emission is beamed or a limited region 
of the ejecta is irradiated, the total irradiated luminos¬ 
ity to the ejecta is smaller than that estimated from 
the observed X-ray luminosity. Note that the ratio 
T'hon/T^rec oc t~ 3 - 88 /t~ 3 oc t~ 088 is decreasing, so that 
the neutral fraction remains unity, / n = 1. 


2.2. Thermalization 

We consider condition Q for the thermalization in the 
ejecta. The optical depth for optical-infrared photons is 
given by 


tir 


kir 


M e j 
47rR ej ' 


(14) 


The dominant opacity at the infrared wavelength is 
bound-bound opacity. In the iron ejecta, the opac¬ 
ity for t he therm al photons is Km ~ 0. 1 cm 2 g -1 
(Kasen, Badnell fe Barnesll2013l:lTanaka &; Hotokezakal 

120131) . For simplicity, we neglect the dependence of the 
opacity kir on the temperature. From condition d4]) and 
Equations 0,0 and m, the upper limit on the ejecta 
velocity is 


— < 2.2 x 10 -2 
c 


ttlR 


0.1 cm 2 g 


2 o—l 


X 


( M ei 

\10~ 2 Mq 


1/2 


1/2 



(15) 


For the r-process eject a,_ (k jr ~ 10 cm 2 g" 1 ; 

Kasen. Badnell fc Barnesl 120131 iTanaka fc Hotokezaka I 
20131 ). the upper limit on the ejecta velocity becomes 
large. 


2.3. Diffusion 

Thermal photons can escape from the ejecta if condi¬ 
tion ([5]) is satisfied. The diffusion timescale bjiff for the 
propagation distance R e j is 


-Rej 

idiff ~ Tir-■ 

c 


(16) 


Then, condition © and Equations (JTJ, ©, (HU) and (fl6l) 
give the lower limit on the velocity, 


— >4.8 x 10~ 4 
c 


KIR 


0.1 cm 2 g 



x 


M, 


ej 


10~ 2 Mq 



(17) 


2.4. Temperature 

Since the upper limit on the temperature was obtained 
from the observations at optical and infrared band s 
(|Tanvir et al. I 120131 iBerger. Fong fe Chornockl 12013l) . 
the temperature of the ejecta has to satisfy condition 
(0 to reproduce the observations Q. The temperature of 
the ejecta is determined by the relation 

4tt R%aT A ~ L m , (18) 

7 Although the lower limit on the temperature is also obtained 
by usi n g the upper lim it on the flux density in 6.7 GHz presented 
by I Fong et al. 1 <120141) , the derived constraint on the ejecta velocity 
is trivial (u e j < c). 
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where a is Stefan-Boltzmann constant. Then, using con¬ 
dition © and Equation (fT8]> . the lower limit on the ejecta 
velocity is 


— >4.1 x 10- 2 
c 


Try 


4000 K 


x 



-2 


Fir 


10 41 erg s _1 


1/2 


(19) 


Since the temperature is determined by the internal en¬ 
ergy density due to the irradiation, the lower limit on the 
velocity m does not depend on the ejecta mass. 


2.5. X-ray and infrared Luminosities 

In Figure HI the X-ray luminosity is approximately an 
order of magnitude larger than the infrared luminosity 
at t ~ 10 6 s. The difference between the X-ray and 
infrared luminosities may come from the collimated X- 
ray emission due to t he geometry of the sup er-Edd ington 
accretion disk (e.g., iTakahashi fc Ohsuea ll2015i l. Note 
that the ejecta velocity is non-relativistic, so that the 
reprocessed infrared emission is isotropic. 

The anisotropy of the ejecta could also reduce the 
infrared luminosity compared to the X-ray luminos¬ 
ity. In the BH-NS me rgers, t he dynamical ejecta are 
highly anisotropic (e.g . iKvutoku. Ioka fc Shibatall2013l 
IKvutoku et al. I 120151 : iKawaguchi et al. 1120151) . If the 
ejecta cover ~10% of the solid angle of the X-ray emis¬ 
sion, the luminosity of the reprocessed emission from the 
ejecta is roughly expected to be Fir ~ O.lFx- 

The anisotropy of the ejecta composition could also 
contribute to the difference between X-ray and infrared 
luminosities. Numerical simulations including the neu¬ 
trino transport showed that the electron fraction around 
the orbital axis is relatively higher than that near the 
equatorial plane (iSekiguchi et al. I 120151: iMartin et al. I 
120151 : iRichers et al. 1 120151) . Then, iron-rich ejecta may 
only concentrate on the orbital axis. In the other re¬ 
gion, r-process elements may be synthesized and hide 
the reprocessed photons with the high opacity. This ef¬ 
fect effectively reduces the emission area of the ejecta 
so that the thermal luminosity becomes lower than the 
X-ray luminosity. 

Conversely, the infrared luminosity may be larger than 
the observed X-ray luminosity because of the limited en¬ 
ergy band of the observations. If a significant energy is 
emitted as unobserved soft X-rays, the infrared reemis¬ 
sion has larger luminosity than the observed X-rays. 

3. RESULTS 

In the left panel of Figure© we show the allowed pa¬ 
rameter space (red region) for the iron ejecta in the X- 
ray-powered model. We use f n = 1, A = 56, T = 4 x 10 3 
K, t = 7 days, «ir = 0.1 cm 2 g _1 and Fir = 10 41 erg s _1 . 
The luminosity, temperature and timescale are consistent 
with the observed values in GRB 13060 3B dTanvir et al. I 
120131: iBerger. Fong fc Chornock 11201.'ll ). Conditions (fT3]l . 
(usd, uni' and (j!9| are plotted as solid lines. This 
panel shows that even if the ejecta are composed of iron, 
the X-ray-powered model can explain the observed in¬ 
frared excess by the reasonable range of the parameters 
(M e j > 4 x 1O~ 2 M 0 and u ej /c ~ 0.04 - 0.1). 


The right panel of Figure [3] shows the case that the 
ejecta is dominated by the r-process elements. The X- 
ray absorption condition requires such as the recombi¬ 
nation rate and photoionization cross section of the re¬ 
process elements. Here, we assume that r-process ele¬ 
ments are partially ionized (not all f-shell electrons are 
ionized) and most X-ray photons are absorbed by the 
ejecta within ~ 10 days. Then, the bound-bound opac¬ 
ity is significant for the optical-infra red ph otons . We 
use ktr ~ 1 0 cm 2 g -1 (fkasen. Badnell fc Barnes 11201,'! 

iTanaka fc Hotokezaka 1120131 ). Other parameters are the 

same as in the iron ejecta case in the left panel of Figure 
[3] The right panel of Figure [3] shows that even if the 
ejecta mass is M e j ~ 10 -3 — 1O~ 2 M 0 , the observed in¬ 
frared excess is reproduced in our X-ray-powered model 
with the ejecta velocity v e -Jc ~ 0.04 — 0.2. 

In the right panel of Figure El we also plot the al¬ 
lowed region for the r-p roce ss nu clear heating model in 
iBerger. Fong fc Chornock 1 (I2013H (blue region). This 
panel shows that the allowed region of the parameter 
space (red region) is much larger than that of the re¬ 
process heating modeQ Since the r-process heating 
model gives the peak l uminosity at the time t ~ Fiiff (e.g., 
iLi fc Paczvhski I f1 99 8l) , the allowed region of parameters 
resides near the line derived from condition (©. The 
X-ray-powered model (red region) has lower ejecta mass 
than the r-process model, so that the contribution of 
the r-process heating is negligible in the X-ray-powered 
model. 

In the both panels of Figure [3] we plot dashed lines 
where the kinetic energy of the ejecta E^ j n equals to the 
internal energy of the ejecta injected by the irradiation 
Flint ~ Tir£ at the time t = 7 days. If the kinetic en¬ 
ergy of the ejecta before the irradiation is smaller than 
the internal energy provided by the irradiation (the grey 
region in Figure©, the ejecta would be accelerated and 
the both energies reach the equipartition, E^ n = E- lnt . 
This means that the ejecta velocity after the irradiation 
is determined by the internal energy deposited by the 
irradiation Flint ~ Fir! and the ejecta mass M e j, 


1 2Fi R t 


M P 


( 20 ) 


as long as the velocity is non-relativistic. We should 
use this velocity for deriving four conditions in the case 
Flkin < Flint. In the r-process ejecta, the allowed 
region with the ejecta mass M e j < lO _3 Af 0 satisfies 
Flkin < Flint- If we consider the ejecta with the mass 
M e j < 1O~ 3 M 0 and the initial velocity v e \ w hich satis¬ 
fies conditions tir > 1 (USD and T < T max (fl9l) . the ejecta 
velocity after the irradiation, Equation exceeds the 
upper limit from condition tir > i ED- Therefore, the 
lower limit on the ejecta mass is M e j ~ 1O~ 3 M 0 in the 
case of the r-process ejecta. On the other hand, for 
the mass range 1O _3 M 0 < M e j <4x 1O _3 M 0 , even 
if the initial ejecta velocity v e j is smaller than the limit 
from T < T max , the irradiation accelerates the veloc¬ 
ity to the allowed region. Note that E[ nt ~ Lir£ is 


8 See IKisaka, Ioka Sz Takamil (120151) for more precise calcula¬ 
tions of the allowed parameters in the r-process heating model. 

9 In the range shown in Figure[3] most of the regions E^ n = Ei n t 
(dashed line) are excluded by the temperature condition (I19D . 
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Fig. 3.— Allowed parameter space to reproduce the infrared excess by reprocessing the X-ray excess emission in GRB 130603B. Left 
panel shows the iron ejecta case. A reasonable range of the parameters, the ejecta mass M e j > 0.04 Mq and velocity v e j/c> 0.04, satisfies 
our model criteria. The condition > ft i on is satisfied in the displayed area with = 10 keV and fix ~ 10 41 erg s 1 . so that f n = 1 
is a good approximation (Equation l|9|)). The kinetic energy is larger than the internal energy above the grey region [the boundary (the 
dashed line) is described by Equation (1201) 1 . which is the physical situation. Right panel shows the r-process ejecta case. The ejecta mass 
A-/ e j ~ HP 3 — 10“ 2 Mq and velocity u e j /c ~ 0.04 — 0.2 satisfy our model criteria. For comparison, we also plot the parameter region of the 
r-process heating model in lBerger, Fong fc Chornock I (120131 ') as a blue region. 


the lower limit on the injected internal energy because 
in the early stage the larger energy may be injected, 
which gives negligible emi ssion at the time ~ 1 0 day s due 
to the adiabatic coo ling dKisaka. Ioka fe Takami 1120151 : 

iKisaka. Ioka &; Nakamura II2015T) . 


4. ON THE OTHER MACRONOVA CANDIDATES 
4.1. Application to GRB 060614 

The model of X-ray-powered macronovae could also 
applicable to the macronova candidate following GRB 
060614. As is the case in GRB 130603B, GRB 060614 
has the long-lasting (~ 10 6 s) X-ray components whose 
timescale and flux are comparable to those of infrared 
excess components (iMangano et~aL 120071 Yang et al. 
12013 1.Tin et ani20laT The luminosity, temperature and 
timescale are ~ 10 41 erg s _1 , ~ 2200 — 3400 K and ~ 12 
day i n the source rest frame, respectively (jYang et al. I 
120151 iJin et al~l 120151) . Because of relatively long 
timescale (~1.7 times longer than that of the macronova 
in GRB 130603B), large ejecta mass M e j ~ O.lAfc is 
requi r ed in the r-pro cess heating model (jYang et al. I 
120151 I.Tin et al ~1 1201, 5j ). Such ejecta mass may be too 
large f or NS-NS mergers a s sugges ted by th e simulations 
(e.g., iHotokezaka et al.l l2013bD . Thus, lYang et al. I 
( 20151) and Jin et al. I (12015 ) sug g ested a BH-NS merger 
( Kvutoku. Ioka fe Shibatal 20131: iKvutoku et al. 1 120151) 


for the origin of GRB 060614. On the other hand, if 
we apply our X-ray-powered model to the observed ex¬ 
cesses in GRB 060614, the required minimum mass is 
~ O.IMq for the iron ejecta and ~ 2 x 10“ 3 Mq for the 
r-process ejecta, where the thermalization and photon 
escape conditions m and 23 give the scaling M e j oc t 2 
compared to GRB 130603B. The required ejecta velocity 


v e - 3 /c ~ 0.03 — 0.0(23 is the same as GRB 130603B (from 
condition ijTojl ). Therefore, a NS-NS merger can explain 
the observed macronova in GRB 060614 more naturally 
in our X-ray-powered model. 

4.2. Application to GRB 080503 

The X-ray-powered model may also explain the opti¬ 
cal rebrightening at ~ 1 — 5 day after the GRB 080503. 

The long-lasting X-ray (~_1 — 2 day ) has a similar 

flux to the optical one dPerlev et al. I I2009D . Unfor¬ 
tunately, the redsh ift of GRB 080503 is only limited 
to z < 4 dPer gy el al. i j jQOO i. From the observations 
(jPerlev et al. 1 120091) . the optical flux and temperature 
are ~ 10 -15 erg s _1 and ~ 5 x 10 3 — 10 4 K, respectively, 
which are roughly constant between ~ 1 — 5 days in the 
observer frame. To explain this event based on the X-ray- 
powered model, the thermalization condition (fl5l) at ~ 5 
days and the escape condition 23 at ~ 1 day have to 
be satisfied. We find that the small mass M e] ~ 10” 4 Mq 
and the velocity u e j/c ~ 0.01 — 0.04 are required in the 
case of the 7'-process dominated ejecta. The redshift is 
also restricted to z < 0.1, otherwise the kinetic energy of 
the ejecta is smaller than the observed radiated energy, 
which is unphysical. On the other hand, the iron ejecta 
with the mass M e j ~ 3x 10~ 3 — 10 _2 Mq and the velocity 
v e j/c ~ 0.01 — 0.04 can also explain the observed optical 
rebrightening in GRB 080503 within the redshift range 
z < 0.5. The allowed parameter region for GRB 080503 
with the iron ejecta slightly overlaps with that for GRB 
130603B with the r- process ejecta, but not with the iron 
ejecta. If GRB 080503 and GRB 130603B have the same 
ejecta mass and velocity, our results imply that the ejecta 

10 Although the velocity is slightly small compared with the 
escape velocity of the neutron stars (~ 0.2c), the uncertainties of 
the temperature which come from the assumed afterglow model 
may reduce the difference. 
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composition is not isotropic. Namely, some parts of the 
ejecta are dominated by iron while other parts are dom¬ 
inated by r-process elements. This picture is supported 
by the numerical simulations (e.g., |Sekigu chi et al~l 
120151: iMartin et, al. II 20 T 5 I : iRichers et, al. 1120151) . 

5. SUMMARY AND DISCUSSION 

We propose the X-ray-powered macronovae as an al¬ 
ternative interpretation of the observed infrared ex¬ 
cess in the short GRB 130603B. Our model has 
important a dvantages over the r -process heating 
model (e.g., iLi fc Paczyhski I I1998H a nd the previous 
engine-powered models (|K saka. Ioka fo Takami I 120151 
iKisaka. Ioka fc Nakamura! 12015T) . First, the X-ray- 
powered model explains both the X-ray and infrared ex¬ 
cess components with similar fluxes using a single energy 
source. On the other hand, the previous models must 
consider different unrelated sources to explain the excess 
components. Second, the X-ray-powered model allows 
for a much broader parameter region of the ejecta prop¬ 
erties than that of the r-process model (in the right panel 
of Figure ED. In particular, the X-ray-powered model 
explains the macronova with smaller ejecta mass. The 
r-process model must satisfy all the conditions which 
the X-ray-powered model satisfies, in addition to hav¬ 
ing the correct amount of r-process material to gen¬ 
erate the observed infrared luminosity by radioactive 
decay, as well as to having the correct diffusion time 
at the observed peak time t ~ tdiff- Then, the al¬ 
lowed parameter range of the ejecta mass has to be 
much narrow in the r-pr ocess model. In the previous 
engi ne-powered model s (IKisaka. Ioka fa Takami I 120151 
IKisaka. Ioka fc Nakamura! 12015ft . the energy injection 
occurs at 10 2 — 10 4 s, earlier than the macronova emis¬ 
sion. To confine the injected internal energy in the ejecta 
up to the macronova timescale ~ 10 6 s, the ejecta sat¬ 
isfy the condition t < tdiff until ~ 10 6 s. This condition 
always requires the larger ejecta mass than that of the 
X-ray-powered model. 

We also apply the X-ray-powered model to the other 
macronova candidates, GRB 060614 and GRB 080503. 
The r-process ejecta with the mass M e j ~ 2 x 10~ 3 — 
3 x 1O _2 M 0 and the velocity v e -Jc ~ 0.03 — 0.08 can 
explain the optical-infrared excess in GRB_060614. As 
oppose to the r-p rocess heating model (lYang et al. 1120151 : 
l.lin et alH 120151 b the X-ray-powered model allows rela¬ 
tively small ejecta mass expected for the NS-NS mergers 
(e.g., iHotokezaka et al.~ll2013bt h In GRB 080503, the 
allowed parameter range does not overlap with that in 
GRB 130603B expect for the case that GRB 080503 has 
the iron ejecta while GRB 130603B has the r-process 
ejecta. This may imply ejecta that are in part made of 
iron and in part made of r-process elements, as suggested 
by numerical simulations that find anisotropic electron 
fraction in the ejecta of NS-NS mergers. 

In the case that the temperature in the ejecta becomes 
low (< 2000 K), dust may for m in the ejecta depend¬ 
ing on the ejecta composition ({Takami. Nozawa fc Ioka I 
120141 ). The dust enhances the opacity of the ejecta. 
Then, the thermalization occurs for lower mass density of 
the ejecta than FigureE] (condition (lT5l) h On the other 
hand, photons can escape in the ejecta only with lower 
mass density (condition (fl7ll l. As a result, the allowed 
region moves toward lower ejecta mass than Figure [3] 


iMetzeer fc Piro 1 (1201411 suggested that ejecta is fully 
ionized and X-ray photons freely escape after ~ 1 day 
for their magnetar model with the X-ray luminosity 
Lx ~ 10 45 erg s _1 . In our model, the irradiated X-ray 
luminosity is much smaller than their magnetar model 
(Lx ~ 10 45 erg s" 1 ), so that the ejecta are not fully 
ionized. 

As the X-ray excess component, we consider the activ¬ 
ity near the central engine, such as the accretion disk. 
Using the fallback model with the temporal evolution of 
the mass accretion rate M oc f -5 / 3 at t > to and the 
radiative efficiency ~ 0.1, the luminosity from the disk is 
~ 10 41 erg s" 1 (M f /O.O17M 0 )(f o /O.l s) 2 %7 day)- 5 / 3 , 
wher e Mf is the total fallback mass (e.g., IKisaka fc Ioka I 
1201511 . Note that the ratio of the fallback mass Mf to the 
ejecta mass M e j is ~ 1 — 10 in our model, w hich is consis¬ 
tent w it h the numerical simulat i ons (Hotokezaka et, al.~1 
l2013allbl : iNaeakura et al. 1 12014 iKvutoku et al. 1120151b 
If the central engine activity comes from the mag¬ 
netar, the luminosity is ~ 10 41 erg s _1 (Rd/2 x 

10 16 G) 2 (R/0.45 s ) -4 where Rd is the surface dipole mag¬ 
netic field an d the P is th e rotation period of the NS at 
7 days (e.g., lUsov Ifl992j ). Since the X-ray excess com¬ 
ponent roughly follows the temporal index ~ —2, the 
period P has been slowed down at 7 days. 

There are important implications of the X-ray-powered 
macronovae for the search of the electromagnetic coun¬ 
terparts to GW sources. If we can first observe 
the X-ray emission originated from the central engin e 
(iKanner et al. l201 2fc IKisaka. Ioka fc Nakamura! 120151) , 
the thermal emission from the ejecta could be estimated 
from the X-ray emission as the luminosity 

L ~ L x , (21) 

and the temperature 



as long as the ejecta satisfy the conditions (ED — ©. Even 
if the ejecta initially prevent the detection of X-ray pho¬ 
tons from the central engine, the ejecta become optically 
thin to the X-ray tens of days after the merger. The 
luminosity is about ~ 10 4 ° — 10 41 erg s _1 , which could 
be detectable by Swift/XKT with inte gration time 10 4 s 
within 100 Mpc ({Moretti et, al. l[2007t ). 

If the X-ray from the central engine is significantly 
faint, other heating sources would become dominant such 
as the r-process and the shorter timescale activities of the 
central engine (extended and plateau emissions). Thus, 
the X-ray observations are important to distinguish the 
main heating source of the macronovae. 
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